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ABSTRACT 

A scenario is put forth for the formation of supermassive black holes at the 
centers of galaxies. It depends upon the formation of a Mestel disk with a flat 
rotation curve, M <r oc r and £ oc 1/r. Such disks could form from the collapse of 
uniformly rotating, isolated, gaseous clouds, either proto-galactic, galaxy-mass 
damped Lya clouds or the gas that survives galaxy mergers. We propose that 
in either case the disk will be unstable to the Rossby vortex instability (RVI). 
This instability grows from any large, steep pressure gradient in an optically 
thick disk. Such pressure gradients either occur adjacent to compact objects 
or could be triggered by heating from individual supernovae in and around 
the disk. Upon excitation, the RVI transports angular momentum outward, 
accreting nearly all mass within the initiation radius. We have calculated that 
in very thin disks, the non-linear vortices initiated by the RVI can transport 
angular momentum far more efficiently than turbulence. Compared to a 
viscosity-based Shakura-Sunyaev disk, the RVI transports angular momentum 
out to a much larger radius, so more mass is accreted into the central black 
hole. A typical galaxy rotational velocity is v rot = 200 km s -1 , and the critical 
column density, necessary to initiate the RVI is Hccd — 100 g cm -2 . For 
M <r = 27rr 2 E, we have t C cd = v^ ot / (2tiT, C cdG), and the mass accreted becomes 
M BH = ^ oi /(27rS CCD G 2 ) = 3 x 1O 7 M . Both the black hole mass M BH and 
its vf ot dependence are in good agreement with recent observations, because 
Vrot = \/3 cr c , where a c is the velocity dispersion of the bulge at the radius of 
mutual contact. 
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1. INTRODUCTION 

There is mounting evidence that non active galaxies, including our own, harbor central 
supermassive black holes (SMBHs) with masses from 10 6 to 10 9 M (Richstone et al. 1998; 
Tremaine et al. 2002). Furthermore, observations indicate a close correlation between the 
SMBH mass and the host galaxy bulge mass, and an even tighter correlation between 
the SMBH mass and the bulge velocity dispersion, a c . How and when the SMBHs form 
has remained unknown, although there are many possibilities (Rees 1984; Ostriker 2000; 
Adams, Graff, & Richstone 2001; Loeb 1993; Eisenstein & Loeb 1995; Umemura, Loeb, & 
Turner 1993). 

Mestel (1963) developed a model of cloud collapse which preserves specific angular 
momentum. His spherical cloud collapsed to a disk with a nearly flat rotation curve, but 
a uniformly rotating inviscid cylindrical cloud should collapse to a disk with a perfect 
flat rotation curve. This mass distribution, M <r oc r, should extend inward to where 
pressure forces become important, e.g., the mass enclosed within a solar radius r becomes 
M = O.2M . The Helmholtz cooling time at every radius is ~ 10 4 yr, which is short 
compared to the BH formation time of 10 8 yr, so the disk is thin. Turbulence has been 
invoked to form a Shakura-Sunyaev (SS) disk at larger radii. However, for a 10 8 M BH, 
an SS accretion disk would become subject to the self-gravitating instability (SGI) at 
~< 0.01 pc (Begelman, Frank, & Shlosman 1989; Shlosman et al. 1990; Goodman 2003; 
Shakura & Sunyaev 1973, Gammie 2001). At 0.08 pc, the disk would have more mass 
than the BH because of the excessive SS column density required to transport angular 
momentum at a rate corresponding to nearly the Eddington luminosity. Hence, even if 
the SGI turbulence is excited and maintained, the disk mass is still greatly excessive. The 
large-scale vortices of the Rossby vortex instability (RVI) disk are one solution to this 
problem. 

An Eddington limit (1 M yr" 1 ), Thompson opacity SS disk surrounding a central 
mass of Mbh — 1O 8 M results in Hss/r — 10~ 4 (Shakura & Sunyaev 1973; Frank, King, 
& Raine 1985; Goodman 2003; Papaloizou & Pringle 1984). Therefore, with the critical 
column density (CCD) necessary to initiate the RVI at t C cd — 10 3 r 55 ~ 10 pc, the mass 
enclosed becomes 7rr 2 S CCD ~ 1O 8 M . All turbulent viscosity mechanisms, no matter how 
robust, will lead to a disk mass too large for the disk to be consistent with the observed 
mass, growth time, and luminosity of active galactic nuclei. The RVI disk does not rely on 
turbulence but on vortices with transport lengths of / ~ rj 10, which is much greater than 
l~H~ r/10 4 for the SS disk. 
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2. THE COLLAPSE OF A GASEOUS CLOUD AND FORMATION OF A 

MESTEL DISK 

A self-gravitating cloud is parameterized by its mass, M v = (47r/3)pr v 3 , initial (virial) 
radius, r v , initial (virial) velocity dispersion, v v = (M v G/3r v ) 1//2 , and the dimensionless spin 
parameter, 

A= JIE^G^M-^ 2 . (1) 

Here J is the total angular momentum, and \E\ — J M <r (G/r)dM describes the potential 
energy (Peebles 1969). For a lone test particle with tangential velocity vy as it orbits a 
fixed mass, equation (1) gives A = v^/vk = 1, where vk is the Keplerian velocity We are 
interested in the final collapsed state of a rotating cloud. The mass at its outer boundary 
will form the outer boundary of the equilibrium disk. However, simulations calculate the 
mass average of A, A at) , not its value at the outer surface. For our purposes we renormalize 
A to \ eq so that a particle in a Keplerian orbit tracing an object's equator has \ eq = 1. 

We also greatly simplify Mestel's analysis by treating a sphere as a right circular 
cylinder with the sphere's mass, uniform density, and equal total angular momentum. 
The cylindrical approximation introduces only a 3% change in radius, yet the expected 
variation in precollapse cloud shapes is very much greater. The advantage of the cylindrical 
approximation is that in Mestel's model it forms a flat rotation disk, because each cylindrical 
shell at r\ can be treated as a separate test particle that collapses to its own Keplerian 
radius, tk- For a cylinder rotating as a solid body, the specific angular momentum starts 
as j(r,j) oc rf and M(< r,j) oc rf , so j(r;) oc M(< Tj). After collapse, j(r K ) oc (Mr^) 1 / 2 , 
so M oc (Mr k) 1 ^ 2 , and one obtains M oc r K : i.e., a flat rotation curve as seen in the 
calculations of Bullock et al. (2001) and Cen et al. (in preparation). 

Using our equatorial definition of \ eq , a typical gaseous cloud, collapsing with no 
transport of angular momentum, becomes rotationally supported at r ~ ^l q r c ioud — 1 kpc 
starting from a typical cloud with r doud = 300 pc, v v ~ lOOkm/s, and M v ~ 10 12 M Q of dark 
and (~ 15%) baryonic matter. A typical value of \ av for a cloud formed in cosmological 
structure formation is 0.05 (Warren et al. 1992; Steinmetz & Bartelmann 1995; Cen et al. 
(in preparation)). The radius of support is about 1/10 of a typical L* galaxy radius of 10 
kpc, a well-recognized dilema, leaving in doubt this explanation of BH formation based on 
the transition from inviscid to vortex flow at the critical column density. However, we see 
two factors that should reduce this large discrepancy: (1) the conversion from \ eq to \ av 
and (2) the reduction of interior mass as the dark matter and baryonic matter separate 
during the collapse. 

For a sphere, the mass averaged \ av = 0.31 if the equator rotates at the Keplerian 
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velocity. This is because an integration of A combines the moment of inertia of a sphere, 
I = (2/5MR 2 u), the self-energy \E\ = (3/5)M 2 G/R, cu = \E\^ 2 /R, and v sph = uR. Thus 
the disk radius predicted by r above should be larger by (1/0. 31) 2 = 10.4, giving agreement 
with observed disk size because of this correction alone. However, a further correction may 
be made for the decrease of total interior mass as the baryonic matter falls deep inside the 
dark matter Navarro- Frenk- White distribution, p(r)/p crit = S c /((r/r s ) ■ (1 + r/r s ) 2 ). For 
r << r s , the interior mass of dark matter decreases as M dk oc r 2 , but the baryonic Mestel 
disk follows Mb ar y oc r, giving pdark « Pccd- We note that the collapsed radius depends 
sensitively on the initial spin parameter \ eq , which implies that the total amount of gas and 
hence the SMBH masses could have a large dispersion, corresponding to the dispersion in 
A e <?- 

2.1. Making a Gaseous Disk by Galaxy Mergers 

In addition to the collapse of a single gas cloud as described above, in galaxy-galaxy 
mergers the gas fraction of the progenitors should collide and make a uniformly rotating 
gas cloud. When the two clouds collide, they produce a gas pressure corresponding to the 
kinetic energy of collision. This occurs before their subsequent collapse by cooling. The \ eq 
of the now combined cloud will depend sensitively upon whether the collision was prograde 
or retrograde (Van den Bosch et al. 2002), but in either case a merged cloud will rotate as 
uniform body, collapse, and form a disk with nearly flat rotation. 

A typical bulge might have 10 10 M or 10% of the mass of a massive galaxy. However, 
because we expect mergers to take place after initial galaxy formation, the gas fraction 
should be 10 9 M or ~ 10% of the bulge stars. The outer boundary of a typical bulge 
merges with the flat disk at ~ 500 pc where the velocity is Keplerian. The relaxation of a 
uniformly rotating bulge of stars and gas should leave the merged stars in inclined orbits, 
thus forming the visible bulge. The gas will collapse, shock, and cool just as in the collapse 
of the original cloud, forming a sub-disk of 10 9 M with a flat rotation curve within ~ 500 
pc. Barnes & Hernquist (1991) and Mihos & Hernquist (1996) have simulated such a merger 
producing a more compact mass inside ~ 200 pc. 

3. THE RVI IN A MESTEL DISK 

The RVI is a global hydrodynamic instability in thin disks, excited by a radial 
extremum in an entropy-modified version of potential vorticity (Lovelace et al. 1999; 
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Li et al. 2000). Its nonlinear evolution has been studied using global two-dimensional 
hydrodynamic simulations in which the flow of matter through the vortices is well illustrated 
in Fig. 6 of Li et al. (2001). One may rightfully ask why the SGI is not sufficient to 
transport the angular momentum as Gammie (2002) has discussed, but if the RVI works, 
then the strong shocks may disperse slowly growing self-gravity perturbations. In fact, the 
SGI might even be useful, because it could help initiate the RVI. The RVI, once initiated, 
has been shown to efficiently transport angular momentum outward, especially by the 
large-scale vortices that it produces. In order to initially excite the instability, there must 
be a large enough local pressure gradient associated with the potential vorticity extremum. 
The RVI can be initiated provided the pressure gradient condition is met. Also, the column 
density of the disk has to be high enough so that locally generated heat can be trapped 
for at least several orbits before radiative cooling damps the fluid velocity stresses. This 
thickness condition is actually quite generic for many hydrodynamic instabilities in disks 
since they rely on pressure forces to drive fluid motions. The only hurdle is the need for 
positively correlated Reynolds stresses for outward transport of angular momentum. 

The growth rate of the RVI depends on the sound speed, but if the collapse were to 
reach equilibrium at a temperature of ~ 10 K, the sound speed would be c s ~ 10~ 3 t> v . This 
low sound speed would combine with an orbit time of t or ut — 3 x 10 5 yr at 10 pc to yield 
an RVI growth time of ~ 10 9 yr, which is unrealistically long. However, during the initial 
collapse there is both a time and a radius at which the column density is great enough to 
trap a major fraction of the heat of collapse, allowing for rapid growth of the RVI. This 
would ideally lead to a sound speed c s = {^ a diab — l) 1 ^ 2 — v v /3, where 'jadiab = 4/3. A 
disk with E = 100 g cm -2 collapsing at 1 pc radius is likely to be initially supported by 
radiation pressure when collapsing through a disk height H pa O.lr, with a corresponding 
T ~ (pc 2 /a) 1/4 = ( (w v / 10) 2 E/alT) 1/4 = 740 K. At this temperature the opacity is of the 
order of k ~ 1 cm 2 g" 1 , and the radiation pressure is ~ 10 x nkT. Consequently, the 
cooling time becomes t coo i ~ H(E/k)/(c/3) ~ 10 _2 fi _1 , far too fast for the RVI to be 
initiated. However, as the disk cools to 100 K, the height shrinks to H pa O.Olr, the disk 
becomes supported by matter pressure below T = 340 K, and with an opacity of k ~ 1 
cm 2 g -1 , the cooling time increases to r cool ~ /J(S//t)/(c/3))(TC p p/aT 4 ) ~ 0.4r2 1 . If we 
consider somewhat greater thickness at correspondingly smaller r, then the cooling time 
scales as r cool ~ #(E//t)/(c/3))(TC>/aT 4 ) ~ [rjr im f^ x = (E/Eioo) 3 ^ 1 . Thus there 
is a thickness 100/t < I^rvi < 1000k g cm -2 where the RVI should maintain itself with 
an instability growth time oc (r/r 10 o) 3 , which is a fraction of the SMBH formation time. 
Three-dimensional simulations are planned, a laboratory experiment has been proposed 
(Colgate & Buchler 2000), based on evidence that Rossby vortices dominate the transport 
of angular momentum in Earth's atmosphere. 
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3.1. Angular momentum transport and triggering of the RVI 

Li et al. (2001) have studied various initial disk pressure profiles that are unstable to 
the RVI in a disk where £ is uniform. (Since the initiating vortices have radii r v <C r, we 
expect the Mestel disk, with X oc r, to behave similarly.) The RVI can have a growth rate 
of ~ 0.1 — 0.3f2 when the initiating scale length associated with the local pressure gradient 
P s /(dP s /dr) is approximately 3H (Fig. 6, Li et al. 2001). Large pressure gradients produce 
strong nonlinear vortices with strong shocks. In these calculations the shocks reach pressure 
ratios of P s > 3pc 2 with c s = O.lvxep- The pressure contrast of such shocks is large compared 
to that expected in strong turbulence where c tur b < 0.3c s , so P shock / Pturb — 3c 2 /c 2 urb ~ 30. 
Thus we postulate that the RVI generates the necessary pressure contrast in shocks 
such that the large-scale / ~ r/10 transport of angular momentum persists regardless of 
how thin the disk is. The torque produced by the RVI becomes H RV i — ij- 2 • H)(P s /3). 
The ratio of torque produced by the RVI to that produced by turbulent viscosity is 
H RVI /H turb ~ 0.03P s /(pc 2 s )(H/r), where H turb ~ (r 2 • H)(pc 2 s /9)(H/r). Both azimuthal 
pressure gradients P s /3 and p(c s /3) 2 • (H/r) act with a torque arm r over an area Hr, but 
the turbulent stress is smaller because the turbulent eddy scale is limited to H (Pringle 
1981). Thus in the special case of a massive BH's accretion disk where H/r ~ 10~ 4 , the 
RVI torque may be greater than the SS turbulent torque by a factor of ~ 10 5 . 

During the cloud collapse, we expect the SMBH to be seeded by a BH from a giant 
central star. Each new stage of the star's evolution causes a steep pressure gradient at 
the inner boundary of the disk. If the RVI starts on the inside, it must progress to the 
outside where all the mass resides. There should also be a significant supernova rate, 
perhaps one per 100 yr per 1O 8 M or roughly 10 6 supernovae in the formation time of the 
SMBH. A single supernova of 10 51 ergs near the disk is sufficient to heat a local area of 
the disk (~ 0.0l7rr 2 E) and raise the sound speed to c s ^sn — O.lt^, high enough to initiate 
and sustain the RVI. In addition, we expect the RVI to initiate limit cycles as mass builds 
up and the instability condition is exceeded. This is strongly reminiscent of cataclysmic 
variables. 

The disk should return to a stable, more quiescent state when sufficient radiation loss 
causes strong damping, ultimately stabilizing the instability at some fraction of the critical 
column density. This is suggestive of what is observed in the water maser nebula NGC 
4258, which is so accurately interpreted as a cool nebula in Keplerian orbit around a SMBH 
(Claussen, Heiligman, & Lo 1984; Henkel et al. 1984). Its column density is estimated to 
be SjVGC — 10 g cm -2 . We therefore interpret NGC 4258 as the residual disk left over after 
forming the SMBH. Such a disk should be self-gravitating, and indeed the drifts in radio 
frequency lines are interpreted (Bragg, Greenhill, & Moran 2000) in terms of spiral density 
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waves. This supports our assumption that the strong shocks of the RVI should suppress the 
SGI until these shocks terminate. 



4. BH MASS AND VELOCITY DISPERSION DEPENDENCE 



In our simplified model of forming a flat rotation disk, a cloud with an initial radius 

\ eq will collapse to r = \ 2 eq i 



rj, total mass M , and spin parameter \ eq will collapse to r = \ 2 eq ri with a column density 
S = Mo/nrl, or S = M /ir(\l q ri) 2 . Then the total mass M C cd within the critical column 
density T>ccd of such a disk is 



MccD = ^cc M D k q n? (2) 

There should be a large dispersion in SMBH masses since the variance to mean ratio 
of A from simulations is <j\ ~ 0.4 (Warren et al. 1992), and Mbh oc MqA -4 . We note that 
for a cloud with uniform density and uniform rotation, any sub-region has the same A as 
the whole cloud, and so the collapse of sub-regions will follow equation (2). 

The rotation velocity of the disk v ro t is related to the disk column density by 
v 2 ot = 27rGS ro, where G is the gravitational constant. From equation (2), the mass within 
a fixed column density T, C cd of such a disk is 

M(> S CCD ) = 11 v? ot » 3 x 1O 7 M , (3) 

lTlG z h C CD 

if v rot = 200km/s and T>ccd — ^rvi — 100 g cm -2 . On very general grounds, we expect 
v rot — \/S a c , where a c is the velocity dispersion of the bulge. This is because at the radius 
where the spherical bulge meets the flattened disk, r^, the interior mass is identical for both 
the bulge and the disk. The extent to which the Keplerian velocity at this point is greater 
than the flat rotation curve velocity at much larger radius is a measure of the increase 
in central mass during mergers, AM merger = M Fiat(( v b I 'v Fiat) 2 — !)• In many cases this 
fractional increase in velocity near the bulge is small, so the increase in mass is small. Hence, 
the mass inside Rccd becomes M<ccd — 2.7 • 1O 8 M , which agrees well with observations 
by Tremaine et al. (2002), who discovered M BH>obs = 1.35 x 10 8 ((T c /200km/s) 4 - 02±a32 . 
For galaxies that form their BH before mergers, we must use v rot and assume that 
the rotation curve is flat all the way to the Rossby radius of ~10 kpc when the 
galaxy is very young. Ferrarese (2002) has discovered a tight correlation between a c 
and v rot : \ogv rot = (0.84 ± 0.09) log a c + (0.55 ± 0.19). Combined with Tremaine's 
observations, this gives M BH = 1.82 • 10 7 (f rot /200) 4 - 79 . This is close to our model, in which 
M BH = 2.88 x lO 7 (*w/2OOA;m/s) 4 M . 
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5. DISCUSSION AND CONCLUSIONS 

We have proposed a simple model for the formation of a SMBH, which is synchronous 
with the formation of its host galaxy from the gravitational collapse of a massive cloud. We 
rely on two assumptions: (1) the initial local conservation of angular momentum leading to 
a flat rotation disk and (2) the subsequent initiation of the transport of angular momentum 
at a given column density of the disk. Together, these lead to a simple explanation of the 
previously puzzling mass of the SMBH and its mass-a relation. We have related the critical 
column density for the onset of angular momentum transport with our work on a large-scale 
instability in accretion disks, the RVI. Star formation may follow exactly the same scenario 
in molecular clouds as SMBHs in Lja clouds. There is growing evidence that gravitational 
tidal torquing of "cores" in star formation regions of molecular clouds lead to the same 
values of X av ~ 0.05 as in Lja clouds for galaxies. Then the mass of the cores within the 
RVI limit becomes Hccd^^ccd ~ l-^o f° r r core — 0.04 pc and \ av = 0.05. 

Gaseous disks formed at the centers of galaxies are not exactly flat rotation disks even 
after gas collisions following mergers. Furthermore, mergers of binary SMBHs are expected 
to evolve somewhat differently because of differing mass ratios and impact parameters. 
These effects as well as the dispersion in \ av will introduce a finite dispersion in the bulge 
velocity, as the M BH — a c relation indicates. This cosmological explanation of the M BH — a 
relation will be explored in a subsequent paper deriving the SMBH mass distribution from 
the \ av and Press-Schechter distributions. 
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